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ABSTRACT 


Applying a transformation to a non-Gaussian field can enhance the information content of the 
resulting power spectrum, by reducing the correlations between Fourier modes. In the con¬ 
text of weak gravitational lensing, it has been shown that this gain in information content is 
significantly compromised by the presence of shape noise. We apply clipping to mock con¬ 
vergence fields, a technique which is known to be robust in the presence of noise and has been 
successfully applied to galaxy number density fields. When analysed in isolation the result¬ 
ing convergence power spectrum returns degraded constraints on cosmological parameters. 
However substantial gains can be achieved by performing a combined analysis of the power 
spectra derived from both the original and transformed fields. Even in the presence of realistic 
levels of shape noise, we demonstrate that this approach is capable of reducing the area of 
likelihood contours within the — erg plane by more than a factor of three. 


1 INTRODUCTION 

The extraction of useful information from cosmological fields 
(such as galaxy number densities or weak lensing) is often heav¬ 
ily compromised by two limiting factors. Firstly the smallest phys¬ 
ical scales are too challenging to model theoretically and are there¬ 
fore excluded from the analysis. This is particularly problematic 
for more evolved fields where the amplitude of fluctuations are suf¬ 
ficiently large that the particle dynamics behave in a highly non¬ 
linear manner. Secondly, for any non-Gaussian field, the power 
spectrum fails to retain all of the statistical information captured 
by the field’s configuration. The amplitude of the field’s Fourier 
modes are correlated, thereby compromising the information con¬ 
tent of the power spectrum. These correlations can be interpreted as 
originating from the largest dark matter halos, which form spikes 
in the density field. Each spike in real space translates to a uni¬ 
form signal in Fourier space, making a coherent contribution to the 
power spectrum across all wavelengths. Both of these issues, the 
modelling of small scale clustering and the correlation of Fourier 
modes, are predominantly associated with the highest density re¬ 
gions of the field. 

Manipulating the observed field prior to evaluating the power 
spectrum can help alleviate these problems. Suppressing fluctua¬ 
tions within the highest density regions in the field improves the 
performance of perturbation theory, allowing a greater number of 
Fourier modes to be included in a likelihood a nalysis. This has been 
demonstrated with simulated galaxy fields dSimpson et alj|201^ 
and applied to data fro m the Galaxy And Mas s Assembly (GAMA) 
galaxy redshift survey dSimpson et alj 12011) . The suppression of 
non-Gaussian peaks also leads to a reduction in the correlation be¬ 


tween Fourier modes, thereby enhancing the amount of informa¬ 
tion retained by the power spectrum. While previous applications 
of clipping involved three-dimensional fields, in this work we turn 
our attention to weak gravitational lensing which involves a projec¬ 
tion of the cosmological density field onto two dimensions. 


The apparent alignment of galaxies on the sky arises from 
the optical distortion imposed by the intervening distribution of 
inhomogeneous matter. When considering the two-point statistics 
associated with this weak gravitational lensing signal, a promi¬ 
nent degeneracy emerges between two cosmological parameters: 
the amplitude of linear density perturbations erg and the cos¬ 
mic matter density Hm. For example, the cosmic shear corre¬ 
lation function from the Canada-France-Hawaii Telescope Lens¬ 
ing Survey (CFHTLenS) has been found to measure the combi¬ 
nation CTg (f2m/0.27)° ® = 0.79 ± 0.03 within the context of 
a standard flat ACDM cosmology jKilbinger et aLlIlQiAh . There 
are a number of possible routes to improve t his measurement. 
Using either tomography i He^mans_et^ 2013h or the full three- 
dimensional information jKitching et alj2014h could lead to signif¬ 
icant gains. Alternatively, further information can be extracted from 
the field by analys ing the three-point statistics, as demonstrated 

' ^ 


by IFu et alj ll2014 ), or t he abundance of peaks jCardone et ^ 
20151: IWhite etal. l2002l: iBerge et alj 1201 ol: iMaturi et alj l20ld: 


Dietrich & Hartlar l2010h . 


ISeo et al] ( 1201 ih found that by invoking a logarithmic trans¬ 
formation, considerable improvements could be made to t he infor¬ 
mation content of the convergence power spectrum, while lYu et alj 
( l2012h demonstrated that further enhancements were achievable 
by performing a wavelet decomposition of the convergence field. 
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2 Simpson et al. 


Ijoachimi et alj ( 1201 ih present a comprehensive analysis of apply¬ 
ing Box-Cox transformations (among which the logarithmic trans¬ 
forms form a subset) to the convergence field. They also reported 
that significant enhancements to the power spectra were achievable 
under idealised conditions. However each study concluded that re¬ 
alistic levels of shape noise substantially reduce these gains. 

Our work aims to build on these previous studies in two ways. 
First of all we shall focus on the clipping transformation which is 
well suited to noisy fields. Secondly, we will investigate the po¬ 
tential benefits of combining the information from both spectra, 
before and after the transformation, by taking into account their 
cross-covariance. 

In ij^we briefly review how local transformations may be ap¬ 
plied to weak lensing fields. Methods for estimating the covariance 
matrices via numerical simulations are detailed in followed by 
our main results which are presented in Sj4] Our conclusions are 
summarised in SjS] 


2 TRANSFORMING WEAK LENSING EIELDS 

Clipping is a local transformation characterised by the application 
of a saturation value (5o to a scalar field S(x) such that 

Sc(x) = So {S{x) > 5o) 

Sc(x) = S(x) (S(x) ^ 5o) ^ 

yielding the clipped field 5c(a:). The motivation for applying clip¬ 
ping to cosmological fields is twofold: (1) the potential improve¬ 
ment in modelling smaller scale clustering, and (2) the introduc¬ 
tion of new information to the power spectrum. In this work we 
shall focus on the latter. 

This transformation has the unique property of either eliminat¬ 
ing or unperturbing the existing noise field. Consider a noise field 
n, defined as the difference between the original field S(x) and ob¬ 
served field Sn (x) 

n(x) ^ Sn(x) - S(x), (2) 

then the transformed noise field n'(x) is the difference in the trans¬ 
formed fields 

n{x) = f{5n)- f{S), 

= n{x)f'{5 ), 

where /' is the functional derivative of /. For clipping, as specified 
by ([D, f'{S) can only take two values: zero or unity. Therefore, the 
noise contribution is either removed (in regions above the thresh¬ 
old) or left untouched (in regions below the threshold). 

The benefits of decorrelating the Fourier modes within the 
galaxy power spectrum are limited since its amplitude is dictated by 
the uncertain manner in which galaxies trace the underlying matter 
distribution. Conversely, enhancing the weak lensing power spec¬ 
trum would be of greater consequence as it acts as a direct measure 
of the amplitude of matter fluctuations. 

A possible disadvantage of applying clipping to lensing stems 
from the fact that lensing involves a projection of the three di¬ 
mensional density field onto two dimensions. This averaging pro¬ 
cess will reduce the level of non-Gaussianity in the field, and so a 
greater proportion of the field must be subject to clipping in order 
to achieve the desired result. 

While the direct observable from weak lensing surveys are the 
shapes of galaxies, it is the derived convergence field which directly 
relates to the projected mass density. The convergence field k{6) is 


therefore the more natural quantity to work with. It may be thought 
of as a weighted integral of the matter perturbations along the line 
of sight 


^ I <?(x)/x(x)^^dx, (4) 

where x denotes the radial coordinate distance. The lensing effi¬ 
ciency (;(x) is determined by the radial distribution of source galax¬ 
ies Ug (x) and the comoving angular diameter distance /k (x) 

,{x) = I ■ ( 5 ) 

If we decompose the matter field into linear and non-linear com¬ 
ponents, 5 = (5g + where 5g is the underlying Gaussian ran¬ 
dom field and the residual term incorporates all departures from 
Gaussianity, then from Eq. the convergence is separable in the 
same manner, yielding k(0) = k.g{0) -f In the limit that 

the non-Gaussian component (6) only contributes at field values 
larger than the chosen threshold, the application of clipping as de¬ 
fined by Eq. ([D ensures that only the Gaussian component of the 
density field Sg{0) contributes to the resulting convergence power 
spectrum . 

Eor comparison we sh all also assess the performance of th e 
logarithmic transformation dSeo et al.ll201 iLlJoachimi et alj|201 ih 

k{0) = Ko In [1 -I- k{9)/kq] , (6) 

where the parameter kq can be adjusted to compensate for the am¬ 
plitude of the convergence field, serving a similar purp ose to the 
clippin g threshold Sq. We follow the prescription used in ISeo et^ 
( I 2 OI ih and set kq = |Kmin| -b 0.001, where Kmin is the most neg¬ 
ative value of the convergence. This small offset ensures that the 
transformed field remains finite. 


3 METHODOLOGY 


In order to quantify the benefits of clipping the convergence field, 
we perform a Eisher matrix analysis to determine the error forecast 
in the Gm — ag plane for a fiducial survey, and compare the size of 
the expected likelihood contours before and after the convergence 
field is subject to a local transformation. 

It is well known that in circumstances where the posterior is 
not well described by a multivariate Gaus sian, the Eisher matri x 
acts as a poor estimator ( see for example ISellentin et alj ( l2014h : 
ISellentin & SchafeJ (12015 1). However when considering forecasts 
for weak lensing surveys in the Gm — erg plane, then provided 
the likelihood contours are sufficiently small, comparisons between 
Eisher predictions and Markov Chain Monte Carlo (MCMC) out- 
puts are found to be high ly consistent (s ee for example Ei gure 3 of 
iTavlor & Kitchindboid and Eigure 6 of lWolz et al.ll2012h . 


3.1 Simulations 

We make use of mock convergence fields from simulations de¬ 
scribed in lHarnois-Deraps & van Waerbeka ( l2015h . We apply a va¬ 
riety of clipping thresholds to the mock convergence fields asso¬ 
ciated with 70 independent lines of sight, each of which spans 60 
square degrees. Sources are placed at a redshift of 2 = 0.8 and 
the fiducial amplitude of linear density fluctuations is taken to be 
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as = 0.826. These lines of sight serve as our estimator for the co- 
variance of the angular power spectra Pnit)- In addition to these, 
we utilise a set of five simulations which are seeded with an initial 
particle distribution with matching phases, but with small displace¬ 
ments in the values of flm and erg. This allows for a clean extraction 
of the partial derivative of the observables with respect to these pa¬ 
rameters, as required for estimating the Fisher Matrices, without 
being contaminated by cosmic variance. Other cosmological pa¬ 
rameters such as the spectral index ris and global curvature are 
taken to be fixed, as they can be measured with a much greater 
precision from other observations such as the Cosmic Microwave 
Background. 

For each line of sight the convergence field is assigned to a 
300 X 300 grid, such that each grid cell corresponds to approxi¬ 
mately 1.5 areminutes, and the convergence power spectra are eval¬ 
uated in eight equally spaced bins for multipoles £ < 1, 500. Our 
results are largely insensitive to the choice of smoothing length. 
However if very fine grids are used then clipping predominantly 
acts on noise spikes rather than true peaks in the projected density 
field. 

When generating each simulated convergence field, the source 
redshift plane ought to remain fixed. However the ray-tracing takes 
place across an integer number of simulation boxes, and since the 
redshift-distance relation is altered when perturbing Hm, there is 
inevitably a small unwanted shift in the redshift of the sources.We 
compensate for this by rescaling the amplitude of the convergence 
field, which accounts for the change in lensing efficiency (;(x) due 
to the small difference between the desired and simulated distances 
to the source plane. The magnitude of this correction is approxi¬ 
mately 1% and so would not change our results by more than a few 
percent if it were neglected. 


3.2 Shape Noise 

Starting from the true convergence field, we use the inverted Kaiser 
Squires method jKaiser & SQuireslll993h to recover the shear field 
at each cell location. Then noise is superposed by drawing two ran¬ 
dom numbers at each cell location, from a Gaussian distribution 
with mean of zero and standard deviation (Je/x/^, which constitute 
the real and imaginary parts of the mean intrinsic ellipticity of the 
galaxies within that cell. The cell averaged ellipticity de is given 
by y/n where n is the number of source galaxies per cell, and 
we take the rms intrinsic ellipticity = 0.3 throughout this work. 
The noisy shear field is then given as a sum of the true and intrinsic 
shear fields. The noisy convergenc e field is recovered by ap plying 
the Kaiser-Squires transformation jKaiser & Sauireslll993h to the 
noisy shear field. 

The intrinsic shape noise introduces an additive white noise 
contribution which includes a distortion associated with the win¬ 
dow function of the grid 

^ ^ \ m 

’ 2n \sinc(fc,/2)2 sinc(fcj,/2)2/ ’ 

where and ky specify the elements of the wavevector in the x 
and y dimensions, in units of the Nyquist frequency, and the aver¬ 
aging takes place within the annulus associated with the wavevec¬ 
tor bin. Clipping reduces the noise contribution to the power spec¬ 
trum, since the regions of the field lying above the threshold are 
smoothed. This leads to an estimate of the cli pped shot noise contri - 
bution Pen (k), derived from equation (5) of lSimpson et alj j2013h . 
where the leading order term may be expressed as 



Figure 1. The ratio of the convergence power spectrum from the clipped 
convergence field, to that of the original field, for three different values of 
ag. In this example a clipping threshold is chosen such that on the largest 
angular scales the power is reduced by a factor of two, which typically 
provides a good balance between the removal of non-linear power while 
maintaining a high level of signal to noise. 


Pen{k)~il-fefPn{k), ( 8 ) 

where /c is the proportion of the field which lies above the clip¬ 
ping threshold. For the fields considered in this work, the value 
of fc required to reduce the large scale power by a factor of two is 
around 20%. By comparison, the proportion of a three-dimensional 
gala xy field subject to a comparable clipping strength is typically 
2% jSimpson et alJl^oT^ . 


3.3 Covariance Estimation 

Since estimating the covariance matrix is a critical step in this anal¬ 
ysis, and one which can be highly susceptible to the influence of 
noise, we compare two different methods for estimating the co- 
variance matrix. First we use the Ledoit-Wolf shrinkage estimator 
jCedoit & Wolfll20o3l to evaluate the covariance matrix C associ¬ 
ated with our binned power spectra. 

C = rp-b (1 - (5*)S , (9) 

We use publicly available cod^ to estimate the shrinkage target 
F and the shrinkage constant 5*. The form of the shrinkage target 
is based upon the samples being independent and identically dis¬ 
tributed. The sample covariance matrix S is defined as the ensemble 
average over 70 lines of sight 

= {APr{i)AP:{j)), ( 10 ) 

where AP^ (i) represents the deviation in power from the sample 
mean, and m denotes the particular transformation(s) under con¬ 
sideration. The total number of columns in the covariance matrix is 
given by the product of the number of angular bins in Pk{(.) with 
the number of transformations included in the analysis, in order to 
accommodate their cross-covariance. 

The finite number of fields is a source of sampling noise which 


^ http://WWW.econ.uzh.ch/faculty/wolf/publications/covCor.m.zip 
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Figure 2. Left: Constraints in the Om — erg plane using the angular power spectra [I < 1,500) of the original convergence field (dashed), the clipped field 
(thin solid) and the log field (dotted). The thick set of contours relate to the combined analysis of power spectra from the original and clipped fields. These 
constraints correspond to a survey of 10,000 square degrees. Right: The same format as the left panel, but we now include the effects of shape noise, with a 
galaxy number density of 20 arcmin“^ and an rms intrinsic ellipticity (Te = 0.3. 


will inevitably propagate through to the estimation of our likeli¬ 
hood contours. In order to confirm that the uncertainty in our likeli¬ 
hood contours is negligible, we employ a second algorithm to esti¬ 
mate the covariance matrix. Instead of using a shrinkage estimator, 
we simply invert the sample covariance matrix S before applying 
a rescaling factor which accou nts for the noise-induced bias in the 
amplitude dHartlap et alj|200% 


c-i 


n — p — 
n — 1 



( 11 ) 


where p represents the length of the data vector (either 8 or 16), 
and n is the number of samples. When all 70 fields are used to 
construct the sample covariance S, these two algorithms produce 
highly consistent results, the derived likelihood contours differ in 
area by less than 5%. 

A further validation test for convergence of the covariance ma¬ 
trix involves repeating this procedure with a small subset of the 
available samples. This also proves useful in order to identify which 
algorithm produces more consistent results. For example, when 
using only 20 samples the rescaling method defined by Eq. (ED 
systematically overestimates the area of the likelihood contour by 
more than 40%. By contrast, the shrinkage estimator recovers like¬ 
lihood contours which are consistently within 40% of the area de¬ 
rived from the full sample. This improved stability leads us to con¬ 
clude that the shrinkage algorithm likely offers a superior perfor¬ 
mance. It also predicts slightly more conservative gains when com¬ 
bining the clipped and undipped spectra. For these two reasons, the 
shrinkage algorithm defined by Eq. ID is the one we adopt when 
presenting our main results. 


3.4 Constructing the Fisher Matrix 

The precision with which a parameter may be determined is gov¬ 
erned by the rate at which the observable responds to a change in 
the parameter. The key observable to consider is the power spec¬ 
trum of the convergence field Pk(^). This may be expressed as 




dx, (12) 


where the lensing efficiency g{x) Is defined in Eq. ([D- 

Predicting the power spectrum of a clipped convergence field 
poses a significantly more challenging problem. As is apparent 
from Eq. (HD, the convergence field can be thought of as aris¬ 
ing from the weighted sum of intervening lens planes. Applying 
clipping to a convergence field derived from of a single lens plane 
would possess a powe r spectrum in line with the model presented in 
ISimpson et alj ( l2013h . However each of the stacked lens planes is 
at a slightly different stage of gravitational evolution. We therefore 
resort to using numerical simulations in order to determine the ex¬ 
pected behaviour. In order to define a consistent clipping strength, 
independent of smoothing length, we specify the fractional drop in 
amplitude of the large scale power. We perform an iterative proce¬ 
dure on the threshold value So to converge upon the desired clip¬ 
ping strength. Note that this iteration is not performed for each field 
- it is performed once, then this threshold value is fixed for all re¬ 
maining fields within the ensemble. Furthermore, prior to applying 
any transformation, each convergence field has its mean subtracted. 
This reflects the uncertain normalisation of a given convergence 
field derived from the measurement of a shear field, although this 
is less problematic for fields spanning a large area. 

Our fiducial survey spans 10,000 square degrees, and involves 
a galaxy number density of 20 arcmin“^. For simplicity we only 
consider a single redshift bin, rather than using tomography, and 
the mock survey is idealised in that there are no holes or masks. 

The Fisher information matrix, formally the expectation of the 
Hessian matrix of the log likelihood, may be expressed as 


F = D'^C”^D, (13) 

where C is the covariance matrix corresponding to the various 
bins in Pk{(.) and is estimated following the methodology in 1 13.31 
Where a combined analysis of two spectra is performed, the co- 
variance matrix incorporates the cross-covariance between the two 
spectra. The Jacobian matrix D is comprised of elements Dij which 
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Undipped 

Logarithm 

Clip 40% 

Clip 60% 

CUp 80% 

Undipped 

1 

1.4 

3.1 

3.3 

2.9 

Logarithm 


0.73 

2.6 

2.5 

1.8 

Clip 40% 



0.59 

0.87 

1.5 

Clip 60% 




0.59 

1.2 

Clip 80% 





0.72 


Table 1. The improvement in the Figure of Merit as defined in Eq. GD, for vai'ious combinations of power spectra from transformed convergence fields, 
relative to that of the original power spectrum. The four local transformations - three clipping thresholds and the logarithmic transform - all prove detrimental 
when considered alone, yet yield substantial gains when coupled with the power spectrum of the original field. 


hold the differential of the iih multipole bin with respect to the jih 
parameter, 


_ dP4i) 

’ dp 


(14) 


We expect transformed fields will be susceptible to shape noise and 
so D must be re-evaluated for each transf ormation when n oise is in¬ 
troduced, as discussed in further detail in ISeo et^ l l2012h . We take 
the average power derived from twenty different noise realisations 
in order to ensure that the above derivatives are stable. 

There are two distinct analyses we wish to perform in the 
— (Js plane. The first is the constraint which one would achieve 
from a conventional interpretation of the weak lensing power spec¬ 
trum, without clipping. In this case the determination of the differ¬ 
ential matrix D can be performed analytically. This is not the case 
for the power spectra derived from the clipped convergence fields. 
In this case, to determine the differential matrix D we utilise the 
simulations which were seeded with matching phases in the initial 
density fluctuations. Due to the nonlinear nature of structure forma¬ 
tion, the smallest density fluctuations exhibit a stochastic behaviour 
across the different simulations, but on the angular scales of interest 
we find the derivatives to be numerically stable. 


4 RESULTS FROM SIMULATIONS 

Using the simulations outlined in the previous section, we explore 
the consequences of applying local transformations to the conver¬ 
gence field, and construct likelihood contours derived from the re¬ 
sulting convergence power spectra. 


4.1 Power Spectra 

Figure [T] shows the change in the shape of the convergence power 
spectrum for three different cosmologies, after applying a clipping 
threshold chosen such that the largest angular bin in Ce is reduced 
in amplitude by a factor of two. Previous studies of clipping cosmo¬ 
logical fields conclude that reducing the linear power by between 
40 — 60% strikes a good balance by suppressing higher order terms, 
while maintaining a high level of signal to noise. This therefore 
drives our choice of the threshold value, which for the central value 
((Ts = 0.831) corresponds to kq = 0.026. The three different val¬ 
ues of (T8 generate significantly different responses, illustrating how 
clipping can reveal the degree of non-linearity present in the field. 
For more evolved density fields, a greater suppression of power oc¬ 
curs on small angular scales. This response can be understood by 
considering the original field to be a superposition of its linear and 
nonlinear components. The clipped power is then predominantly 
comprised of the clipped linear power spectrum and the clipped 


nonlinear power dSimpson et al.ll2oT^ . However owing to the pos¬ 
itive skewness of the nonlinear pdf, the nonlinear power suffers a 
much stronger loss of power than the Gaussian pdf. Thus, fields 
which are more evolved and host a stronger nonlinear component 
suffer a greater loss of power on small scales. This behavior has 
previously been obse rved for the case of th e three-dimensional sim¬ 
ulated density fields dSimpson et alj2013h . Applying stronger clip¬ 
ping with a lower threshold value leads to a more linear power spec¬ 
trum with less correlated Fourier modes. This improvement satu¬ 
rates once the amplitude of the linear power spectrum is reduced 
by approximately a factor of two, so there is little to be gained by 
adopting a lower threshold. 


4.2 Parameter Constraints 

The left hand panel of Figure illustrates the 68 and 95% joint 
likelihood contours derived from the Fisher matrix of the conver¬ 
gence power spectrum {£ < 1500). The dashed lines are derived 
from the standard convergence field, while the dotted lines are de¬ 
rived from the log transform as defined by Eq. ill}. The outer set 
of solid contours make use of the clipped convergence field as de¬ 
fined by Eq. (D, with a threshold selected such that the large scale 
power is reduced to 50% of its original value. These contours are 
somewhat elongated relative to the original spectrum, but also ex¬ 
hibit a small clockwise rotation of their degeneracy direction. The 
inner set of solid contours result from combining both the original 
and clipped power spectra. The minor axes of the combined con¬ 
tours are not noticeably reduced from those of the original spec¬ 
trum, owing to the strong cross-covariance between the amplitudes 
of the two spectra. However the major axes shrink considerably, 
suggestive of a powerful complementarity between the information 
held by the two power spectra. 

Correlations between data points are often perceived to have 
a detrimental impact on the dataset’s overall information content. 
But there are many circumstances in which this is not the case. To 
see why, consider the simple case of two data points, to which we 
wish to fit a straight line. The only two free parameters are the gra¬ 
dient and the amplitude. If the errors on these two data points are 
perfectly correlated, the constraint on the amplitude weakens by a 
factor of x/2 relative to the uncorrelated case. However the gradi¬ 
ents uncertainty vanishes - it is now determined to an arbitrarily 
high precision. 

This toy model is intimately linked to the case we consider 
here. For weak lensing contours in the Dm — os plane, the size 
of the minor axis relates to the uncertainty in the amplitude of the 
weak lensing power spectrum. This is often quoted as a measure¬ 
ment of the parameter combination asHm- The fact that the clipped 
and undipped spectra are highly correlated means that when a joint 
analysis is performed, any improvement in the measurement of this 
amplitude is negligible. Hence the minor axis does not shrink for 
the case of the ‘Combined’ contours in Figure 2. Meanwhile the ex- 
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tent of the contour’s major axis is governed by how well the shape 
of the power spectrum is measured. More evolved fields with higher 
values of era will possess more power on smaller angular scales, 
and this helps lift the degeneracy within the flm — era plane. When 
considered separately, the precision with which the shape of each 
spectra can be determined is ultimately limited by cosmic variance. 
It might be that the large scale power is low by chance, causing a 
spurious upward tilt in the lensing power spectrum, thereby mim¬ 
icking a higher value of aa. However this deception can be revealed 
with a joint analysis of both spectra. A higher value of era leads to 
a stronger decrement of small scale power in the clipped spectrum 
relative to the undipped. This relative change in power is largely 
insensitive to the particular cosmic variance realisation. The con¬ 
tour’s major axis therefore shrinks further than would have been 
possible if the two spectra were uncorrelated. This is explored fur¬ 
ther in il4.3l and Figure[3below. 

Note that this is not the only cosmological example where cor¬ 
relations prove to be beneficial. Performing an analysis of redshift 
space distortions with multiple tracers allow cosmic variance fluc¬ 
tuations to be effectively nulled, permitting an enhanced measure¬ 
ment of the growth of structure. _ 

_ These results are consistent with the findings of ijoachimi et alj 

( I 2 OI ih who noted that upon performing a Box-Cox transformation 
the likelihood contours demonstrate a small clockwise rotation of 
the degeneracy direction, and become considerably elongated. This 
elongation of the contours for transformed fields occurs because 
information regarding the shape of the spectra (as induced by non¬ 
linear structure associated with higher values of as) is erased. In¬ 
deed the full extent of this elongation effect is likely lost in our case 
due to the limitations of the Fisher formalism. Nonetheless the ad¬ 
dition of the original information from the undipped field truncates 
the longer contours. The combined contours therefore offer a faith¬ 
ful representation of the true l ikelihood. _ 

Our findings differ from ISeo et alj (1201 2h in that even in the 
noise-free case, we do not observe significant gains from apply¬ 
ing a logarithmic transform alone. However there are a number of 
differences in our methodologies. They use multiple redshift bins 
and explore a broader parameter space - for example we do not 
marginalise over the dark energy parameters wo and Wa- Further¬ 
more the efficacy of the log transform is strongly governed by the 
choice of pixel size. 

The right hand panel of Figure [2] takes the same format as the 
left panel, but here we include the impact of shape noise. There is 
only a modest degradation of the undipped contours in this case. 
The degeneracy direction of the transformed fields experiences a 
slightly reduced rotation, but there is still sufficient complementar¬ 
ity to yield substantial gains when performing a combined analysis. 
The logarithmic contours also show little degradation, which may 
appear surprising given the information loss found by ISeo et al] 
( I 2 OI ih . However the impact of shape noise is highly sensitive to 
the choice of smoothing length used to define the field, particularly 
since the amplitude of the convergence fundamentally alters the 
transformation defined in ([^, and here we are using a considerably 
coarser pixel size. 

The extent to which the benefits of applying clipping and log¬ 
arithmic transformations are degraded by the presence of shape 
noise also depends on the specific choice of transformation cho¬ 
sen. Higher clipping thresholds lead to a higher power spectrum 
and thus a lower susceptibility to shape noise. Yet they are also less 
efficient at recovering information to the power spectrum in the low 
noise limit. 

In order to quantify these gains we can readily compute the 



Figure 3. Simulated constraints for a 200 square degree convergence field, 
before (red dashed contours) and after (blue solid contours) applying a log¬ 
arithmic transformation (95% C.L.). The inner black solid contour repre¬ 
sents a combined analysis of both .spectra. The consequence of ignoring the 
cross-covariance between the original and transformed spectra is illustrated 
by the dotted black contour. Shape noise is not included in the convergence 
field. The black mai'ker signifies the fiducial cosmology of the simulation. 


change in the figure of merit such as those previous ly defined in 
lAlbrecht et al.l ll2009h and lSimpson & PeacockI ( 1201 Oh 


FoM = 


y^FnnFo- 


— F^ 

^ Qa ’ 


(15) 


where the notation Fee denotes the elements of the Fisher matrix 
associated with the parameter 9. This quantity is inversely propor¬ 
tional to the area enclosed by the likelihood contours. The abso¬ 
lute value is not of great interest, since it will strongly depend 
on the configuration and area of the survey in question, so in Ta¬ 
ble [T] we quote in each case the fractional change in the Figure of 
Merit relative to the original power spectrum analysis. We show re¬ 
sults for three different clipping thresholds, corresponding to those 
which reduce the large scale power to 40, 60 and 80% of its orig¬ 
inal value. In addition we include the original convergence field, 
and the log transform as defined by Eq. (|6j. The diagonal elements 
demonstrate that for each transformation under consideration, the 
areas enclosed by the ellipses have enlarged. The off-diagonal ele¬ 
ments in Table[T]quantify the significant gains which can be made 
by performing a combined analysis, particularly when combining 
the original spectrum with that from a clipped field, where fhe gain 
is more fhan a factor of three, despite the presence of shape noise. 
We have also evaluated intermediate clipping strengths, finding that 
the gain derived from combining two clipped power spectra deteri¬ 
orates rapidly as the difference between their clipping strengths is 

reduced. _ 

Thus while the finding s of previous related studies dSeo et alj 
I 2 OIIL Ijoachimi et al.ll2oTI ) are entirely consistent with the results 
presented here, our conclusions differ markedly in that we find sub¬ 
stantial gains can be achieved from a local transformation of the 
convergence field even in the presence of realistic levels of shape 
noise. The source of this difference stems from performing a com¬ 
bined analysis with the original power spectrum, which holds com- 
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plementary information, rather than considering each power spec¬ 
trum in isolation. 

4.3 Full likelihood model 

The posterior in the parameters flm and cts as produced by con¬ 
temporary cosmic shear surveys is not well approximated by a 
multivariate Gaussian. In order to assess whether our conclusions 
are affected by the limitations of the Fisher matrix formalism, we 
consider the case of a ful l mock likelihood ana lysis following the 
methodology presented in ijoachimi et al.l ( 1201 ih . In this work ana¬ 
lytical models for power spectra of convergence fields transformed 
via Box-Cox transformations, which include the logarithm, were 
derived using a perturbative expansion and recalibration of higher- 
order contributions via the third to fifth m oments of the c onver- 
gence measured from N-body simulations. Ijoachimi et alj ( 1201 ih 
used a suite of 100 convergence fields at a single source redshift 
^ = 1 simulated with the SUNGLASS pipeline dKiessling et alj 
I2OI ih for signal and covariance estimation via Eq. GD- We adopt 
these measurements for convergence field without shape noise and 
smoothed with a Gaussian of width 1.5 arcmin. Moreover, we 
rescale the covariance to correspond to a 200 deg^ survey. This 
toy model setup is deliberately chosen to yield wide contours with 
a strong non-linear degeneracy, which is least well represented by 
a Gaussian approximation to the posterior. We use the value for the 
free parameter in the logari thmic transformation, equivalent to kq 
in Eq. (5), as determined by Ijoachimi et al.l 1 I2OI ih in the Box-Cox 
optimisation. 

The resulting parameter constraints for a power spectrum 
analysis of the original convergence fields, of the log-transformed 
fields, and their combined analysis are shown in Eigure[3l We have 
used 10 angular frequency bins logarithmically spaced in the range 
150 to 1500. The contours for the log-transformed fields are slightly 
biased high, which is caused by the difficulty of analytically mod¬ 
elling the transformed power spectr um at high £ where hi gher-order 
correlations become important (see ljoachimi et al]l201ll for a more 
detailed discussion). However, this does not affect our conclusion 
that there is excellent qualitative agreement between the full likeli¬ 
hood analysis and the Fisher forecasts. The dotted contours repre¬ 
sent a joint analysis in which the cross-covariance has been set to 
zero, as would arise if the two spectra were taken from fields in dif¬ 
ferent parts of the sky. This greatly underestimates the achievable 
constraints, demonstrating the importance of the cross-correlations 
in suppressing noise contributions. 

This independent analysis offers striking confirmation of the 
considerable gains which are achievable by accounting for the 
strong cross-covariance which exists between the power spectra of 
a locally transformed field and its original form. In terms of the 
areas enclosed by the contours in Figure [3] the combined analysis 
offers a gain in performance of more than a factor of five compared 
to the transformed field alone. 


5 CONCLUSIONS 

In this work we have demonstrated how applying a transformation 
such as clipping to a convergence field can not only introduce new 
information to the resulting power spectrum, but that this new infor¬ 
mation helps lift the degeneracy in the — os plane. When con¬ 
sidering moderately large angular scales i < 1500, parameter con¬ 
straints can be improved by more than a factor of two, even when 
realistic levels of shape noise are included. The loss of performance 


reported in previous studies for a single transformed power spec¬ 
trum is largely attributed to the elongation of the contours along 
the degeneracy direction, but this elongation is readily truncated by 
including the original power spectrum which exhibits a different 
degeneracy direction. 

A second major advantage of performing a local transfor¬ 
mation such as clipping is that it may also permit smaller an¬ 
gular scales to be studied, since the resulting power spectrum is 
more straightforward to model. We may therefore have signifi¬ 
cantly underestimated the gains achievable when clipping conver¬ 
gence fields. The use of tomography - dividing the source galaxies 
into discrete redshift bins - in conjunction with clipping will further 
enhance the level of precision attainable, but we leave this topic for 
the subject of future work. A suite of N-body cosmological simu¬ 
lations would be required to adequately model the clipped conver¬ 
gence power spectra. Many simulations will in any case be required 
for modelling the conventional power spectrum, or alternatively 
rescaling methods could prove effective for exploring a greater vol¬ 
ume o f parameter space d Angulo & Whitell20ld : lMead & PeacockI 
|2014) . 

As found when applying clipping to large scale structure, the 
choice of threshold is fairly flexible. Provided the threshold is suffi¬ 
ciently strong to reduce the large scale amplitude of the power spec¬ 
trum by over 30%, there is little change in performance. The opti¬ 
mal choice ultimately rests upon the level of noise within the partic¬ 
ular field under investigation. For the convergence fields considered 
in this work we find that in order for clipping to reduce the large 
scale power by a factor of two, a threshold value of kq = 0.026 
is required, corresponding to approximately 20% of the field being 
clipped. 

The greatest obstacle to applying a local transformation to 
weak lensing data is that the observable field is not a direct tracer 
of the matter density, so we must convert the shear field to a con¬ 
vergence field in order to find a suitable proxy for the projected 
matter density. The conversion process is particularly challenging 
in the presence of noise and masks, but significant pr ogress is being 
made ( IVan Waerbeke et^l2013l : IChang et alj|2015| ). We therefore 
envision that a combined analysis of clipped and undipped lensing 
observations will help extract more robust and precise information 
from future and forthcoming surveys. 
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